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' This paper is devoted to the negative flows of the AKNS hierarchy. The main result of this work is 

the functional representation of the extended AKNS hierarchy, composed of both positive (classical) 
and negative flows. We derive a finite set of functional equations, constructed by means of the Miwa's 
shifts, which contains all equations of the hierarchy. Using the obtained functional representation 
we convert the nonlocal equations of the negative subhierarchy into local systems of higher order, 
derive the generating function of the conservation laws and the N-dark-soliton solutions for the 
extended AKNS hierarchy. As an additional result we obtain the functional representation of the 
Landau-Lifshitz hierarchy. 
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1. Introduction. 



This paper is devoted to the negative flows of the AKNS hierarchy. The most well-known 
^ ■ physical models described by this class of equations are the sine-Gordon and self-induced 

transparency models (see, e.g., books [H|9] and references therein). Mathematically, the 
simplest way to describe the subject of this work is to use the zero-curvature representation 
(ZCR), the approach which is the base of the inverse scattering transform (1ST). The AKNS 
hierarchy is an infinite set of equations coming as compatibility conditions for the system 
of linear equations consisting of the Zakharov-Shabat scattering problem [18^119] . 



¥. = U(A)¥, = * A ) (1.1) 

and a linear problem describing the evolution, 

* t = V(A)tf (1.2) 

where V(A) is a matrix polynomial in A. The negative flows that form the negative AKNS 
hierarchy correspond to the case when V(A) is a polynomial in inverse powers of A. 



l 
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Alternatively, negative (sub)hierarchies can be denned in terms of the recursion opera- 
tor [51 lllirnZl ll5j: if the hierarchy of higher-order symmetries (equations of the hierarchy) 
is generated by powers of the recursion operator, then the negative (sub)hierarchy can 
be constructed using its negative powers [15] . The recursion operator technique is widely 
exploited in the so-called "structural" approach [11] to integrable systems when one studies 
such questions as the Poisson structures, multi-Hamiltonian structures, mastersymmetries, 
i?-matrices, etc. 

A characteristic feature of the negative flows is their nonlocality. Among large number 
of known negative hierarchies there are only two, to our knowledge, examples of local ones: 
the Ablowitz-Ladik hierarchy, whose simlest positive and negative equations can be written 

as 



and the Manna-Neveu generalization [TO] of the Hunter-Saxton equation [6 J that was dis- 
cussed in [3]. Looking for the simplest "negative" matrices V(A), V(A) oc A -1 , one arrives at 
nonlocal equations. This nonlocality can be, in principle, eliminated by going to the equa- 
tions that are not of the evolution type (the most well-known example is the sine-Gordon 
equation) or by introducing additional variables (as in equations that appear in the theory 
of the self-induced transparency). However, in any case the standard 1ST scheme is to be 
modified and probably because of this fact the negative AKNS hierarchy is not studied as 
comprehensively as the positive (classical) one. 

The purpose of this paper, where we present some generalizations of the results obtained 
in [2117118]. is twofold. First, we want to describe the extended AKNS hierarchy composed of 
both positive and negative flows. To do this we will use an approach that can be viewed as 
an alternative to the "traditional" ZCR or the method of [2J. The main result of this work is 
the functional representation of the positive and negative AKNS subhierarchies (Sees. [3][5]). 
We derive a finite set of functional equations, constructed by means of the Miwa's shifts, 
which contains all equations of the extended hierarchy, that can be recovered by the power 
series expansion. This approach has some advantages over the standard 1ST because in its 
framework one can avoid introducing the "intermediate" objects of the inverse scattering 
method like Jost functions or scattering data and formulate results explicitly in terms of 
solutions (as, e.g., in the case of the generating function for the conservation laws presented 
in Sec. 

The second goal of this paper is to obtain some explicit solutions, the dark soliton 
ones, of the extended AKNS system. Here we use another advandage of the functional 
representation which in many situations facilitates the calculations, especially when one 
deals with the whole hierarchy (AKNS hierarchy in our case) instead of one of its equations 
(say, the nonlinear Schrbdinger equation (NLSE) in our case). As is shown in Sec.[7J we can 
enhance the classical results of [TO] using simple algebraic proceeding. 

The obtained results have some interesting byproducts. In Sec. El we derive the func- 
tional reprresentation of the Landau-Lifshitz hierarchy (LLH), that is known to be gauge 
equivalent to the AKNS hierarchy, and demonstrate that its negative part is symmetric to 
the positive one, which means that we have another example of local negative subhierarchy. 
Also we show that among the equations of the extended AKNS hierarchy one can find some 




±i 



k = constant 



(1.3) 
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models that were not associated with the AKNS system, such as, e.g., the Wadati-Konno- 
Ichikawa-like equations [16] and Zakharov (1 + 2)-dimensional NLSE [17J. 

2. Holonomy representation of the extended AKNS hierarchy. 

Instead of the ZCR that is based on presenting the equations as the compatibility conditions 
for the linear systems 

A* = U } *, £»-0»., A *_1.2.„. (2.1) 

we will be dealing with the linear systems constructed by the Miwa's shift operators that 
are applied to functions of a doubly infinite number of arguments, 

Q (t, t) = Q (*i,*2> — >*ij*2 5 •••) = Q (*i>*fc)j,fc=i I 2,... ' ( 2 - 2 ) 

and are defined by 

E € Q(t,t) = Q(t + i[f],t), (2.3a) 

l,Q(t,t) = Q(t,t + i[7/]) (2.3b) 

or 

%Q(ij,ffc)i,fc=i,2,... = + it 3 /j,tk)j,k=i,2,...> ( 2 - 4a ) 

1,7^(^,4)^=1,2,... = Qfo,** + iv k /k)j,k=i,2,-- ( 2 - 4b ) 

Thus, our goal is to study the compatibility conditions of the systems of the following type: 

E 5 $ = (2.5) 

and 

1^ = Lfa)¥, (2.6) 

where L and L are 2 x 2-matrices, that are given by 

[E ?1 L(6)]L(ei) = [%L(6)]Lfe), (2.7a) 
jE^Lfa)] !_(£) = [E„L(0]L(»7), (2.7b) 
[\L(l2)] Lfal) = [E n2 L(r/!)] L(r&). (2.7c) 

Using the standard strategy of the ZCR, that consists in introducing an auxiliary parameter, 
£, and looking for the matrices L and L with simplest dependence on £, one can come to the 
matrices that lead to the AKNS hierarchy. It turns out that the matrices L and L should 
be linear functions of £ _1 and £. Indeed, expanding, for example, (|2.5p in the power series 
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in £ one obtains an infinite set of equations of the following structure: 

[d j + V{d^ 1 ,...,d 1 )]t> = L j y, j = l,2,... (2.8) 

where d~ = d/dtj. Expressing recursively the derivatives dj-i, 82 one can rewrite the 
last equation as 

=[i j -v(di,...,d 1 )]^, j = l,2,... (2.9) 

If l_i is linear in £ _1 , then, after replacing the powers of d\ by powers of l_i one arrives at 

d n * = V j 9, J = 1,2, ... (2.10) 

where Uj is a jth order polynomial in £ _1 . In such a way, the standard structure of the ZCR 
of the AKNS hierarchy is reproduced, with £ _1 playing the role of the spectral parameter. 
Omitting the details of the calculations we present here the "minimal" solution of (|2.7p 

and 

CW = 1 -TT^(^ (Mr) (2 - llb) 

where 1 is the unit matrix and the functions Q, R, q and r are subjected to some constraints 
that will be discussed below. 

In the following sections we study the systems of equations (the subhierarchies of the 
extended AKNS hierarchy) that appear as the result of (|2.7|) combined with (|2.11|) . 

3. Positive AKNS subhierarchy. 

The equations that follow from (|2.7ap are 



(6 - 6) %%Q = A (6,6) (6 %Q - £2 %Q) 
(£1 - £2) R = A (£1,6) (£1 - £2 



(3.1) 



with 

A (6, 6) = l + £i£ 2 (%E 6 Q)i?. (3.2) 

By simple algebra one can verify that (|3.ip ensure vanishing of all components of the 
matrix equation (|2.7ap . The functional representation ()3.ip of the positive AKNS hierarchy 
has been derived in |13pi4] (see also [I]). It can be simplified by various limiting procedures. 
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For example, sending £2 to zero one arrives at the system that was used in |13j : 

(Q-Q + i^dxQ - i 2 Q 2 R = 

{ . - (3-3) 

[R-R- i£diR - £ 2 QR 2 = 

where the notation dj = d/dtj and 

Q = E 5 Q (3.4) 
is used. Another reduction can be made by introducing the operator <9(£) by 



Noting that 



one can rewrite (13.11) as 



9(0 ( 3 - 5 ) 

J'=l 



lim -J— fE^El 1 - l) / = iT^CO / (3.6) 



-id(Z)Q = Q-Sl(OQ 
id(£)R = R- R 



(3.7) 



where 



= _ (3.8) 

and 

.R = E^i?. (3.9) 

Expanding these functional equations in the power series in £ one obtains an infinite set 
of differential equations, the first non-trivial of which are the NLSE, 



id 2 Q + d u Q + 2Q 2 R = 
-id 2 R + d u R + 2QR 2 = 



2 _ n ( 3 - 10 ) 



and the complex mKdV equation, 

d 3 Q + d ll iQ + 6QRd 1 Q = 



d 3 R + d lu R + 6QRd 1 R = 

(here stand for d 2 /dtjdt^, etc) that are the first equations of the AKNS hierarchy. 

The positive AKNS sub hierarchy, to repeat, is the classical AKNS hierarchy, that has 
been introduced in the early 1970's and which is one of the best-studied integrable systems. 
That is why we do not discuss Eq. (|3.ip here in details. The only thing that we need 
to illustrate some features of negative AKNS equations is the generating function for the 
conservation laws. 
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Constants of motion. 

The AKNS hierarchy, as an an integrable system, possesses an infinite number of constants 
of motion that can be represented in the form 

lm(t 2 ,t 3 , ...) = I m (ti,t 2 ,t 3 ,...)dti, m = 0,l,... (3.12) 

J reg 

(the symbol "reg" indicates that one has to regularize, if necessary, the integrands by 
adding some constants depending on the boundary conditions that ensure the existence of 
the integrals). As has been shown in [13], the generating function for I m , 

oo 

HO = E ( 3 - 13 ) 

m=0 

has a very simple form when rewritten in terms of Miwa's shifts: 

I(() = (E C Q)R. (3.14) 

Indeed, Eqs. (|3.ip imply that 

d(om=9iJ(c,o (3.i5) 

where 

J{ ^ =(m HmM^ (3 ' 16) 

which leads to djl m = for any j and m. It should be noted that paper [T3] is devoted to the 
NLSE and not to the whole AKNS hierarchy, so a reader can find there only the particular 
case of Eqs. (|3.15p and (|3.16p . However, in order to not deviate from the main topic of this 
paper, we present them here "as is", leaving the proof for the a separate publication. 

4. Mixed AKNS subhierarchy. 

The subhierarchy that is discussed in this section is closely related to the equations that 
usually appear in the works devoted to the negative flows of the AKNS hierarchy. The 
results presented below can be viewed as a generalization of the ones obtained in [TIE] - 
Substituting matrices ()2.11ap and (|2.11bp into (|2.7bp one arrives at 

(r]h(i])q = Q-Q, 
1 r]h(t])r = R — R 



and 

f n — fl _ PnK\ ( n -L 

(4.2) 



q = (1 - (qR) (q + iQ 
r = (1 - iQr) (r + £R) 



where 



m = j^- (4.3) 
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with the shortcuts that are used throughout this section: 

/ = E ? /, f = E v f. (4.4) 

Thus, we have four equations for four functions Q, R, q and r. Since we are discussing the 
AKNS hierarchy, our current task is to eliminate the last two and to obtain a closed system 
of equations for Q and R. 

Equations (|4,ip taken for two negative shifts, and E % , can be rewritten in terms of 
the operator <9(t/), 

d(ri)=Y jV l d l = -i V lim (e^E" 1 -!) (4.5) 



(here d\ = d/dt{) as 



where 



and 



irj d(rj) Q = ui (rj) q, (4.6a) 
-ir/~ 1 d(r]) R = ui (rj) r, (4.6b) 



/ = E„'/. (4. 



Equations ()4.6ap and (|4.6bp can be generalized by means of (|4.ip and ()4.2|) . 

(1 - to) irt^Bir,) Q = u( v )(l- CQf) (q + £Q) , (4.9a) 

- (1 - to) "T^fa) R = w (r?) (l - (r + £R) (4.9b) 

with 

/ = Ef V- (4.10) 
Repeating this trick for two positive shifts, E^ and E^ 2 , one arrives at 

iC'diO = {l-qr- 2^qR) Q, (4.11a) 
1 + qr 1 + qr 

-iC^iO ^— = (l-qr- 2£Qr) R. (4.11b) 
1 + qr 1 + qr V / 

These equations suffice for achieving our goal of eliminating q and r. The resulting system 
can be written as 

f = (1 - to) 0(Od(v)Q + 2i£P(T,)d(£)Q + i[to~ 2£a(£)] d( V )Q + 2f(3( v )Q, 

{ = (1 - to) d(0d(v)R ~ 2itP(r l )d(0R -i[to- 2ta(0] B(j,)R + 2^( V )R, ' "" ~ ' 
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where the functions a and f3 are defined by 

«">=firf (413b) 

and are related by 

d(r,)a(0 = d(OP(v)- (4-14) 

The last identity, that can be verified directly, is very important from the viewpoint of the 
conservation laws of the hierarchy. 

Expanding Eqs. (|4.12p in the power series in £ and rj one arrives at 

f = (d k+ idi + i - d k di + 2iPi +1 d k - 2ia k di +1 ) Q, 

{ _ _ k,l = 1,2,... (4.15) 

[ = {dk+idi+i - d k di - 2i/3i+id k + 2ia k d l+ i) R, 

where a k and are the coefficient of the Taylor's expansion of the functions and /3(f?), 

oo oo 

a (o=j2z k <xk, m = J2 r i l Pi- (4.i6) 

k=l 1=1 

As in the case of positive subhierarchy, this two-parametric system can be simplified in 
several ways. First, Eqs. (|4.12p in the £ — > limit can be represented in the form 

(0 = d 1 d(r ] )Q + tr ] d(r ] )Q + 2(3( V )Q, 

< _ _ 4.17 

[ = did(ri)R - ir]d(r])R + 2^{-q)R, 

where P(rj) can be viewed as an additional dependent variable related to Q and R by 

<9i/3(7?) = B(r])QR (4.18) 

which is the limiting form of (|4.14p . Here one can see the nonlocality (explicit expression for 
f3(rj) invokes d^ 1 operator) that was observed in all works devoted to the negative AKNS 
flows. Equations ()4.17j) can be bilinearized by introducing the tau-functions r, a and p by 

p(rj) = d 1 B(rj)\nT (4.19) 

and 

Q = -, R=£ (4.20) 
r r 

In new terms Eqs. (|4.17p and ()4.18p become 

f = (-Di + iff) D(rj) c • t, 

\ (4-21) 
[0 = (D 1 - i V ) D{r]) p ■ t 

and 

D u t ■ t = 2pa, (4.22) 
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where D\, T)\\ and D(rj) are the Hirota's bilinear operators, 

Dj u ■ v = (dju) v — u (djv) 



D 



jk.. 



DjD^... etc., and 



(4.23) 



(4.24) 



fc=l 



with 



Dk u ■ v = (df.u) v — u (d^v) . (4-25) 
Returning from the power series to "individual" flows one can rewrite (|4.17|) and (|4.21|) 



as 



and 



' = (d 1 Bi +1 + idi + 2(3 l+1 ) Q, 
= [b\B l+l - iBi + 2A+i) R, 
dxft = BiQR, 

f = (DtD l+1 + %Di) a ■ r, 
[0=(D 1 D l+l -iDi)p-T, 



1 = 1,2,... 



(4.26) 



I = 1,2, 



(4.27) 



together with 1^22]) . 

Equations (|4T26|) and (|Q7]) were discussed in 018], thus Eqs. (jlT?) and (jOT|) can be 
viewed as their compact form while Eqs. (|4.12p and ()4. 15|) as their generalization. 

The simplest (and hence most representative) equation of the mixed AKNS subhierarchy, 
discussed in this section, can be written as 



'0 =Q xy + 2PQ, 
< = R xy + 2PR, 

where x = t\, y = i\ and P = /3\. The next one (Eq. (|4.26[) with I = 1) 

[ -iQt = Qxy + 2PQ, 
iRt = R X y + 2PR, 

{ P, = (QR) y 



(4.28) 



(4.29) 



(x = t%, y = t2, t = t\ and P = /?2) is nothing but the (1 + 2)-dimensional NLSE introduced 
by Zakharov [17] . 



Constants of motion. 

In Sec. [3l we have presented the generating function for the constants of motion, /(£) = 
(E^Q) R, for the positive (classical) AKNS hierarchy. It turns out that equation similar to 
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(|3.15p holds for the negative flows as well: 

9(77)I(C) = &LJ(f7,C)> 

where 



'' W (T7)[l-C(E ( g)(E„ 1 r) 



1-Cr? 



(4.30) 



(4.31) 



This means that the quantities I m , as is expected, are the constants of both positive and 
negative AKNS subhier archies. Again, the proof of (|4.30p and (|4.3ip will be published 
elsewhere. 



5. Negative AKNS subhierarchy. 

Equations that follow from the commutativity condition (|2.7cp can be written as 

f Vi 1 {K - 1) HmWrnQ = V2 1 {®m ~ l) H m )E Vl q, 

\ r/f 1 (E m - 1) h( V2 )r = ^ (E„ 2 - l) h( m )r, ' 

where ^(77) is defined by (|4.3p . This is a closed system for the functions q and r which is 
closely related to the LLH that is discussed in Sec. El However, as in Sec. [H we use the 
words "negative AKNS (sub)hierarchy" bearing in mind another system, the one for the 
functions Q and R that can be written in terms of Miwa's shifts and differential operators 
with respect to negative "times", E^ and dj = d/dij. 

There are several ways to eliminate q and r together with from Eqs. (|2.7b|) and (|2.7c|) . 
The shortest one can be described as follows. Passing from E^ to d{rj) one can rewrite (|5.ip 
as 

J id (77) UJ q = Ku(i])q - j(r])u q, 
[ —id (77) UJ$r = Kuj{rj)r — 7(77)0^0^ 

where oj(rj) is defined in (|4.7|) . q = E^q, f = E^ V, 

7(r ] )= 1 —^ = -m, (5-3) 
1 + qr rj 

k = 7(0) = (1 — qr)/(l + qr) and loq = w(0) = 1/(1 + qr). Combining the above formulae 
with Eqs. (|4.ip one can obtain the system 

{0 = d\d(rj)Q + iri^nd^Q - i^^diQ, 
(5.4) 
= 8id(r])R - irT x Kd{jl)R + il{ri)d x R. 

To finish the derivation of the negative AKNS equations one has to express 7(77) in terms 
of Q and R. This can be easily achieved by means of Eqs. (|4.ip : 

7 2 (t?) = 1 - 4T7" 2 [9(t7)Q] [d( V )R] . (5.5) 

So we have a closed system of Eqs. (|5.4p and (|5.5p that can be viewed as the functional 
representation of the negative AKNS subhierarchy. Expanding (|5.4p and (|5.5p in the power 
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series in 77 one arrives at an infinite number of partial differential equations describing the 
negative flows. The simplest one can be written as 

f = (id 2 + ^B^BA Q, 

I (5.6) 

[0= (id 2 - K^iK- 1 ^) R 

with 



, 1-4(B 1 Q)(B 1 R). (5.7) 
The second negative AKNS equation can be written as 

f = (2B 3 - Bm - 3iK 2 B 2 K- 1 B 1 ) Q, 

S - - (5-8) 

[ = (2d 3 - Bm + 3iK 2 d 2 n- 1 di) R. 

One can see that the negative AKNS equations have the same nonlinearity structure as the 
Wadati-Konno-Ichikawa equations [16], that have been derived as a generalization of the 
positive AKNS hierarchy. 

6. Landau-Lifshitz hierarchy. 

In this section we discuss the "pure negative" equations, the ones stemming from (|2.7cp . 
from a standpoint different from that of Sec. To this end it is convenient to introduce 
the matrix 

*">=«m(V-i+*.). <«> 

where q = K v q and r = V. In principle, one can express S in terms of Q and R using ()5.ip . 

S(n)-( 7(7?) - 2 ^ 1B ^) R \ (62) 
S(7?) " {^B( V )Q - 7 (ry) )' (6 ' 2) 

However, we will not use this relationship below, restricting ourselves to the consequences 
of (|2.7cp that can be formulated in the terms of the matrix (|6.ip . By straightforward algebra 
one can show that equations (|5.ip lead to the following equations: 



f]\ 1 B(m) s (%) = % ld (m) s im) (6-3) 

and 

2m / 3 l S(»7) = [S(77),S], (6.4) 

where 

S = S(0). (6.5) 

Using ()6.3p one can rewrite (|6.4p as 

2iB( V )S=[S( V ),S] (6.6) 

and calculate S(ry) in terms of S: 

S(?7) = -iS 5(r))S + X(rj) S. (6.7) 
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Here the function A (rj) should be determined from the condition 

S 2 (r?) = 1 (6.8) 

which leads to 

\ 2 ( V ) = -±tr[d( V )S] 2 . (6.9) 

Thus, we have a closed system of Eqs. (|6.6p . (|6.7p and (|6.9p for the matrix S. 

Expanding these equations in the power series in rj, one can obtain an infinite set of 
equations 

= i&S + i[S",S], (6.10) 

o = a3S + s / " + ^((s') 2 s)', (6.ii) 

where the symbol ' is used to denote the derivative with respect to t±: 

S' = 8iS. (6.12) 

The above equations are the simplest equations of the LLH. In other words, we have shown 
that Eq. (|2,7cp lead to the LLH and derived the functional representation of the latter. 

The fact that the LLH is closely related to the AKNS hierarchy is not new, it is known 
since the works of Zakharov and Takhtadzhyan [20]. However, we would like to note that 
the Landau-Lifshitz equation that was mentioned in [8], whose results are generalized in 
this section, and the Landau-Lifshitz equation that appear in |20] are not the same: the 
last one belongs to the positive subhierarchy, while the former describes the negative flows. 
This indicates that the symmetry between the positive and negative flows of the AKNS 
hierarchy, which is not visible in terms of Q and R, becomes apparent at the level of the 
LLH. 



7. Dark solitons of the extended AKNS hierarchy. 

In this section we would like to present the dark-soliton solutions of the extended (describing 
both positive and negative flows) AKNS hierarchy. We will not derive them from scratch 
but will use the classical results (for the positive subhierarchy) and extend them to cover 
the negative flows. The A^-soliton solutions for the NLSE were obtained in the beginning of 
the seventies by Zakharov and Shabat who developed in |18pi9| the corresponding version 
of the 1ST. Since all equations of the AKNS hierarchy, considered from the viewpoint of 
the inverse scattering approach, are based on the same scattering problem, their solutions 
(dark solitons in our case) possess the same structure that the ones derived in [19]. Thus, 
to solve any equation of the hierarchy one can utilize a big part of the results of [19J. The 
only thing that one has to do is to establish some relations between parameters of the 
solutions (the so-called "dispersion laws") which are different for different equations of the 
hierarchy. These considerations suggest the following procedure: we look for the solutions 
whose structure is similar to the classical dark solitons and then find, using some simple 
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algebraic calculations, the conditions that convert them into solutions of all equations of 
the hierarchy (both positive and negative). 

The main building blocks for the dark-soliton solutions of the AKNS hierarchy are N xN 
matrices A that satisfy the "almost rank-one" condition 

LA - AL _1 = \£){a\ (7.1) 

where L is a constant diagonal matrix, \l) is constant ^-component column, \£) = 
(£i, ...,£n) T and (a\ is ^-component row depending on the coordinates describing the AKNS 
flows, (o(t,t) I = (oi (t, t) , oat (t, t)), and matrices are defined by 

H C = (C1-L)(C1-L- 1 )- 1 , (7.2) 

where 1 is the N x N unit matrix. 

The remarkable property of the above matrices, that will be repeatedly used below, is 
that the determinants 

w(A) = det|l + A| (7.3) 

satisfy the Fay-like identity 

(£ - V) u Zv + (v-() ^ ^vC + (C - Ur, W(£ = 0, (7.4) 

where 

ui = uj (A) , ujq = uj (AH^) , uj£ n = ijj (AH^H^) . (7.5) 

One can find an elementary proof of this identity in Appendix A| It is shown below that 
upon representing the Miwa's shifts as multiplication by combinations of matrices it is 
possible to derive from (|7.4p the identities similar to equations (|3.ip . (|4.ip and (|5.ip that 
we want to solve. 



7.1. Solution of the equations of the positive subhierarchy. 

First let us study the positive subhierarchy. The key feature is to assume that the depen- 
dence on the positive "times" t\ , t 2 , ... is governed by 

E f A = AH Q H 1 , (7.6) 

where the function a = a(£), a(0) = 0, is specified below. Writing down Eq. (|7.4p with 
£ = a\, Tj = ot2 and Q = and the matrix A being replaced with AHq 1 one arrives at 

(ai - a 2 ) w(A) [EiE 2 w(B)] = a-i [Eiw(B)] pE 2 w(A)] - a 2 [Eiw(A)] [E 2 w(B)], (7.7) 

where B = A Ho and 

E fc = E ffc , a k = a{i k ) (k = 1, 2). (7.8) 
In a similar way, Eq. (|7.4|) with £ = 1/ai, r/ = a 2 , C = and A — > AHq 1 leads to 

(1 - aia 2 ) w(A) [E^ 1 E 2 w(A)] = [Ej" 1 a;(A)] [E 2 w(A)j - a x a 2 [E7 ^(C)] [E 2 w(B)], (7.9) 
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where C = AH . Rewriting (|7.7p and (|7.9p in terms of functions Q and R denned by 



Q = U 



;(B) 



i? = y 



c(C) 



w(A)' " ' w(A) ; 
where J7 and V are two auxiliary functions one can obtain 



E12Q 



(E 12 U) pEiw(A)][E 2 a;(A)] 



and 



ai-o 2 £j(A)[Ei 2 w(A)] 
Eiw(A)] [E 2 w(A)] 1 



Ml %Q)--^-(E 2 Q) 



Eit/ 



1 



(Ei 2 t/)y 



E 2 C7 



(E12 Q) R 



w(A)pEi 2 w(A)] l-aio-2 
It is easy to see that these equations become (|3.ip if the following conditions hold: 

(a - 6) « fc E12 u 



(ai - a 2 ) (1 - aia 2 ) E fc 
I aia 2 = (E12 C/) V. 

The simplest way to satisfy these equations is to take 



(fc = l,2) 



E|C7 



f(0 



and 



which reduces (|7.13|> to 



and 



-i/V = constant 



(6 - 6) / (6) / (6) = [1 - (6) / (6)] [&/ (6) - 6/ (6)] • 

The last equation can be transformed, in £ 2 — > limit, into "ordinary" one, 

a2£ a / 2 (0 + (/*£-!)/(£) + 1 = 



(7.10) 



(7.11) 



(7.12) 



(7.13) 



(7.14) 



(7.15) 



(7.16) 



(7.17) 



(7.18) 



with an arbitrary constant /*. Solution of this quadratic equation that satisfies /(0) = 1 
determines the dependence of a on £. 

In a similar way one can show that R defined in (|7,10p satisfies the second of equations 
(pTTj) . Thus, definitions (I7TT0]) together with (fTTHO) . ([7TT5D and ([7TT8D provide iV-dark-soliton 
solutions for the functional equations (|3.ip describing the classical AKNS hierarchy. Below 
one can find a more detailed version of these formulae written down for the physically 
relevant case R = —Q. 
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7.2. Solution of the equations of the negative subhierarchy. 

Assuming 

E„ A = A Hfl 1 H 



A) n P 



(7.19) 



where f3 = j3{rf) and /?o = /?(0) one can obtain from (|7.4|) with ^ = /3, 77 = /?o, C = an d the 

(7.20) 



shift A -> AH^ 1 



(/3 - A))w(A r ) [E„a;(B«)] = /3w(A) [E^(B)] - A)w(B) [E^(A)] . 
In a similar way, Eq. ([73D with £ = P, 77 = l//3 , C = and A ->■ AH^ 1 leads to 

(1 - A,0)w(A) [E^(A)] = u>(A r ) [1^(A«)] - /3 /3^(C) [I^(B*)] . 
Here the matrices A, B and C are the ones defined above while 

B 9 = AH^, C" = AH^, 



A '= AH r/ib. a- = ah 1/a> . 



In terms of the functions 



r = v 



oj(A t 



these equations can be represented in the form 
w(A) [I„u;(A)] 



W ( A -) [E^(A')] E„ g 



and 



(1 - A>/3)w(A) [E„ w(A)] = w(A p ) [!„ w(A«)] 

which leads to (|4.ip after imposing the conditions 

/3 /3 = - (I» 

[U(P-/3 )(l-l3oP)=Vl3a®v u - 
These restrictions can be resolved as follows: 

E V U E v u 



1 ^ (w \ 

1 - 7= , [E v q) r 



(E v u) 



v 



U 



u 



and 



u v 

- = - = s*, 



(7.21) 



(7.22) 



(7.23) 



(7.24) 



(7.25) 



(7.26) 



(7.27) 
(7.28) 



/3 = Pog(v), Po = -uv, 
where g* is an arbitrary constant and g(rj) is the solution of the quadratic equation 

[g( V ) - 1] [1 - $g( V )] = g^girf) (7.29) 

satisfying g(0) = 1. 
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This completes settling the problem of finding iV-dark-soliton solutions of the extended 
AKNS hierarchy because one can show by straightforward algebra that the functions q and 
r presented in this section satisfy Eqs. f)5. If) as well. 

7.3. R = —Q case. 

This section is devoted to the situation that appears in the physical applications of the 
NLSE (and hence, of the whole AKNS hierarchy): 



fl(t,t) = -Q(t,t), (7.30) 

where overbar stands for the complex conjugation. In this case the background solutions U 
and V can be represented in the form 

U = Q* exp [icp (t,t)], (7.31a) 

V = R*exp[-i<p(t,i)] (7.31b) 

with the constants Q* and R* being related by 

R* = -Ql- (7.32) 

Similar formulae can be written for q and r, 

u = g# exp [i(p (t,t)], (7.33a) 

v = r* exp [— i(p (t,t)] (7.33b) 

with 

r* = — qZ. (7.34) 
The phase ip is determined by the equations 



exp[i(E e -lH =/(0, 
exp [i (E v - l) ip] = g(r}). 



(7.35) 



Presenting ip as 



<P (t, t) = cpo + ^ (<^fc*fe + ftktk) (7.36) 
fc=i 

one can obtain for the generating function for the coefficients cp^ and (pk the following 
expressions: 

oo 

J2<PkH k /k = -lnf(0, (7.37a) 

fc=l 

oo 

$^W7*/fc = -]n<7(77) (7.37b) 
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or, after applying the £,d/d£ and rjd/drj operators and using Eqs. (|7.18p and (|7,29p for /(£) 
and g{rj), 



ti i-\Q*\ 2 ep(o 



k [l-g(v)] l~\q*\ g(v) 



k=l 



i - \q*\ 2 g 2 (v) 



(7.38a) 



(7.38b) 



with arbitrary real <pq. Then, it is easy to check that to ensure the necessary properties of 
the matrices A one has to choose 



L = diag(e ie ")r i .=i,...,7v- 
In this case the matrices describing the t- and t-evolution are unitary, 

H 1 H Q(?) =diag(e^«)) 

V / n=l,...,I 

and 



,N 



A ~ & y e / n =l...N 



with 



</>„(£)= 2 arg l-|Q*|e/(CK 



(7.39) 
(7.40) 
(7.41) 
(7.42a) 



MO = 2arg 1 - \q*\g(rj)e- Wn ). (7.42b) 
Now one can establish the dependence of the matrix A on the variables tf. and 



A(t,t) = A* diag 



n=l,...,N 



where A* is a constant matrix and 



V n (t, t) — V nQ + { v nktk + ^nkh 
k=l 



with arbitrary real v n Q and 



k=l 

oo 



fc=l 



or 



fc=l 

oo 



(7.43) 
(7.44) 

(7.45a) 
(7.45b) 

(7.46a) 
(7.46b) 
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Upon noting that the determinants (|7.3p are invariant under transformations A — > M _1 AM 
one can eliminate, without loss of generality, the constants £ m (£ m — > 1) by redefining 
the functions a n (t,t) (a n (t,t) — > a n (t,t)£ n ) thus arriving at the final expressions for the 
N-dark-soliton solutions of the extended AKNS hierarchy: 

Q(t,t) = |Q*|e^)^M, (7.47) 

1\q (t,t) 

R(t,t) = -\Q,\e-^^§. (7.48) 

Here, the determinants A e (e = 0, ±1) are given by 

exp (y n (t,t) + 2ieO n ) 



A £ (t,t) = det 



sm I 



(7.49) 



m,n=l N 



while the phase of and the real constants a n (0, 0) are respectively absorbed into ipo and 
Appendix A. Proof of (I7.4|) . 

In this appendix we present some identities for the matrices defined in Sec. [7] which provide 
a proof of the Fay's identity (|7.4p . Consider the matrix A satisfying 

LA-AM = |£)(a| (A.l) 

with arbitrary diagonal matrices L and M together with the matrices defined by 

H C = (C-L)(C-M)- 1 , (A.2) 

where (£ — L) stands for (£1 — L) etc. It follows from (jA.ip that 

(C- M)(1 + H C A)(C- M)- 1 = l + A-\£){b ( \ (A.3) 



with 



which leads to 



with 



and hence to 



(b c \ = (a\(( - M)- 1 (A.4) 
det |1 + AH C | = det |1 + A| • det |1 - |e)(6 ? | | (A.5) 
|e) = (l + A)- 1 ^) (A.6) 



^ = l-(6 c |e). (A.7) 

A little bit more cumbersome calculations lead to the following "two-point" analogue of 
COD: 

(£ - M)(r? - M) (1 + H S H„A) (£ - M)" 1 ^ - M)" 1 (A.8) 

= l + A+k±\£)(b ri \ + 1 p±\£)(bz\ (A.9) 
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and 



det 



1 + \ui + | u 2 ){v 2 \ 



Here, the rows (vi^\ and the columns | U\ >2 ) are defined by 

(vi\ = (b v \, (v 2 \ = (k\ 

and 



^1) 



1 



rj -£ 



1 



with 



|cc) = (l + A)- 1 (C-L)K). 
Rewriting the determinant in the right-hand side of (jA.lOj) as 



l + («i|iti) (yi\u 2 ) 
{v 2 \ui) l + {v 2 \u 2 ) 



and calculating the scalar products, 



(v 2 \ ui) 
(vi\u 2 ) 



with 



and 



£-77 

¥>C = <fc-cl c C> 



1 + {v\\ U\) = — H 



w 77 - £ 
l + (, 2 |n 2 ) = ^ + ^ 



lo f - 77 



(here Eq. (|A.7|) was used) one arrives at 



which leads to "separation of variables" , 



^ _fr ] _ 



where 



(A.10) 

(A.ll) 
(A.12) 

(A.13) 
(A.14) 

(A.15) 
(A.16) 

(A.17) 

(A.18) 
(A.19) 

(A.20) 
(A.21) 
(A.22) 



Upon adding three copies of (|A.21|) for (£,77), (77, £) and (CO one can obtain the identity 
(|7.4p that we want to prove. 
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